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Abstract 

We report a measurement of the branching fractions in the J/^vr+vr" decay based on a 

140 fb~^ data sample collected at the T(45) energy with the Belle detector at the KEKB e"'"e~ 
collider. Charged pion pairs are found to arise mainly from and /2 mesons; an upper limit 
for the non-resonant contribution is also set. The following branching fractions are obtained: 
B{B^ -> J/V'P°) =(2.8 ± 0.3 ± 0.3) x lO'^, B{B^ ^/V'^) < 1-5 x 10-5(at 90% C.L.) and 
B{B^ J/V'(7r+7r-)non-res.) < 1 X lO'^ (at 90% C.L.). 

PACS numbers: 
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The decay 5° ^ J/i)p^ which is governed by the h ccd transition, can exhibit a 
CP- violating asymmetry. Since the tree diagram of this transition has the same weak phase 
as 6 — i> CCS, a measurement of indirect CP violation in this decay provides an alternative 
estimate of sin20i. In contrast to the b ccs case, however, both tree and penguin 
amplitudes contribute to the b ccd transition in the same order of sin^c- Therefore, 
if penguin or other contributions are substantial, a precision measurement of the time- 
dependent CP asymmetry in 6 — >■ ccd may reveal the values that differ from those for 
b — >• CCS. Thus, B decays induced by the b — > ccd transition, such as —>■ J/ipp^, play an 
important role in probing non-tree diagram contributions. 

Since the p° meson has a large width, it is necessary to study all the decays of the neutral 
B meson that result in a J/ilJ7r~^7r~ final state. In general, there are two possible contributions 
to the B^ — > J/tp7i~^7r~ mode. One is resonant, when a n~^7r~ pair arises from some resonant 
state. The 5° J/ipp^ decay mode is the largest in this class, although it is important 
to search for other resonant contributions. The other contribution is non-resonant, where 
a neutral B meson decays directly into a. J/tp and 71+ vr" pair. For this process, the CP 
eigenvalue is unknown; therefore, the evaluation of this contribution is important to control 
uncertainties when measuring CP violation in B^ J/ipp^ decays. 

The first attempt to measure the branching fraction of B^ J/ipp^ was made by the 
CLEO collaboration who set an upper limit 13{B^ J/i^P^) < 2.5 x 10^'^ at 90% C.L. using 
2.39 fb~^ of data P]. Recently the BaBar collaboration reported their measurement of the 
branching fractions B{B^ J/t/'tt+tt-) = (4.6 ± 0.7 ± 0.6) x 10"^ and B{B^ J/i^P^) = 
(1.6 ± 0.6 ± 0.4) X 10"^ based on a data sample of 52 fb~^ In this paper, we report 
a measurement of the branching fractions of the neutral B meson decays resulting in the 
J/ijj7r~^7r~ final state via B^ J/ipp^ and B^ — >■ J/ipf2 as well as an upper limit on the 
non-resonant 5° — >■ J/ijj'K~^iT~ decay mode. These measurements are based on a 140 fb~^ 
data sample, which contains 152 million BB pairs, collected with the Belle detector |^ 
at the KEKB asymmetric-energy e~^e~ (3.5 on 8 GeV) collider |^ operating at the T(45') 
resonance. 

The Belle detector is a large- solid- angle magnetic spectrometer. Closest to the interac- 
tion point is a three-layer silicon vertex detector (SVD), followed by a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like ar- 
rangement of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(TZ) crystals (ECL). These subdetectors are located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An iron fiux-return located outside of 
the coil is instrumented to detect K^^ mesons and to identify muons (KLM). The detector 
is described in detail elsewhere Q. 

Hadronic events are selected if they satisfy the following criteria: at least three recon- 
structed charged tracks; a total reconstructed ECL energy in the center of mass (cms) frame 
in the range between O.ly^ and 0.8^/s, where ^/s is the total cms energy; an average ECL 
cluster energy below 1 GeV; at least one ECL shower in the region —0.7 < cos 6* < 0.9 in the 
laboratory frame; a total visible energy, which is the sum of charged track momenta and to- 
tal ECL energy, exceeding 0.2^/s, and a reconstructed primary vertex that is consistent with 
the known interaction point. After imposing these requirements, the efficiency for selecting 
i?-meson pairs that include a, J/ if) meson is estimated by Monte Carlo (MC) simulation to 
be 99%. To suppress continuum events, we require the event shape variable R2 to be less 
than 0.5, where R2 is the ratio of the second to the zeroth Fox- Wolfram moment 0. 

The B candidate selection is carried out using two observables in the rest frame of the 
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T(4S') (cms): the beam-energy constrained mass Mbc = \J El^^^ — and the energy 

difference /S.E = Y,Ei — Ebeam, where Ebeam = -\A/2 is the cms beam energy, and Pi and Ei 
are the cms three-momenta and energies of the B meson decay product candidates. 

Candidate J ftp mesons are reconstructed via their decay into oppositely charged lepton 
pairs {e~^e~ or fi~^fj,~). Leptons are selected from charged tracks satisfying a cut \dz\ < 5 cm, 
where dz is the track's closest approach to the interaction point along the beam direction. For 
electron identification, the ratio between the charged track's momentum and the associated 
shower energy {E/p) is the most powerful discriminant. Other information, including the 
CDC specific ionization, the distance between the ECL shower and the extrapolated track, 
and the shower shape is also used. Muons are identified by requiring an association between 
KLM hits and an extrapolated track. Both lepton tracks must be positively identified as 
such. In the e~^e~ mode, ECL clusters that are within 50 mrad of the track's initial momen- 
tum vector are included in the calculation of the invariant mass (Mee(7)), in order to include 
photons radiated by electrons/positrons. The invariant masses of e~^e~{'y) and ii^ ii~ combi- 
nations are required to fall in the ranges —0.15 GeV/c^ < (Mee(T,) —Mj/^) < +0.036 GeV/c^ 
and -0.06 GeV/c^ < (M^^ -Mj/,/,) < +0.036 GeV/c^ respectively. Here Mj/^ denotes the 
world average of the J/ip mass [71]. We perform a vertex fit to the lepton pair candidate. 
Then a mass constrained fit is applied to improve the AE and Mbc resolutions of the selected 
B meson candidates. 

Information from the ACC, TOP and CDC is combined into a likelihood ratio for 
kaon/pion separation; we then impose a cut to reject kaons. The efficiency for pions is 
more than 85% while the kaon fake rate is 10%. Charged pion candidates are also rejected 
if they are positively identified as leptons. Using pion candidates with \dz\ < 5 cm, n'^ir' 
pairs are formed. We apply a vertex fit to tt+tt^ pairs and require the distance between 
the reconstructed vertices of the J/tp and the tt+tt^ pair to be less than 3 mm in order to 
suppress B^ —>■ J/ipKg, Kg tt+tt" events and backgrounds due to accidentally formed 
pion pairs. 

J/tp candidates and tt+tt^ pairs are combined to select B candidates. The distribution of 
the candidates in AE and Mbc is shown in Fig. ^ The signal box to select candidate events 
is defined as 5.270 GeV/c^ < Mbc < 5.290 GeV/c^ and -0.04 GeV < AE < 0.04 GeV. 
Using a MC data set of B^ — > J/ipn^n^ decays distributed uniformly in phase space, the 
detection efficiency is estimated to be 27.8%, where a difference of the pion identification 
in efficiency between MC and data is taken into account: the MC efficiency is 2.8% higher 
per pion track than in data. After applying all selection criteria, there are 537 candidates 
remaining in the signal box. 

A binned maximum likelihood fit is performed to the distribution of the two pion invariant 
mass (M^+^-) of the selected events to determine various contributions to 5° J/tp7i^n~ 
events. Five types of events are considered: (i) 5° — > J/ipp^; (ii) B^ J/ "^12', (iii) 
B^ —>■ J/ipTT^n' (non-resonant signal); (iv) B^ — > J/ipKglKg —>■ 7r"^7r~); (v) background. 
A probability density function (PDF) is constructed for each of these five cases. The B^ — »• 
J/ijjKg mode is not considered as a signal while determining the branching fractions for the 
decay modes contributing to B^ —>■ J/tpTT^n^. 

The PDF used to model the B^ J/ipp^ mode is a relativistic P-wave Breit-Wigner 
function 0]: 

F(M ) - M...-r(M...-)P^^e«+^ 
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FIG. 1: (a) Distribution of AE versus M^c for — > J/ip7r~^7r candidates in data, (b) Projection 
in AE for 5.270 GeV/c^ < Mbc < 5.290 GeV/c^ and (c) Projection in Mbc for -0.04 GeV < A^; < 
0.04 GeV. The arrows in (b) and (c) indicate the boundaries of the signal box. 



where 

- r„(^)-H5|i^),i±|M). (2) 

Here, g(M^-i-^-) is the pion momentum in the di-pion rest frame, with go = Qi^p)] P is the 
J/ip momentum in the rest frame; / = 1 is the p° meson's spin; Mp = 775.8 MeV/c^ and 
To = 146.4 MeV/c^. Also, Les is the effective orbital angular momentum between the J/ij) 
and the which can take any value between and 2, and R is the Blatt-Weisskopf barrier- 
factor radius [9|]. The fit is performed with L^s and R equal to 1 and 0.5 fm, respectively. 

The PDF for the 5° — > J/ipf2 mode is a relativistic D-wave Breit-Wigner function that 
is obtained by replacing the appropriate parameters in the PDF for J/ipp^: I = 2, 

= 1.285 GeV/c2 and Fq = 184.3 MeV/c^. 

The PDF used to model the 5° J/ipKg mode is a single Gaussian function with 
the mass and width fixed to the values obtained by fitting a Kg — >■ n~^7r~ invariant mass 
distribution in 5° J/ipKg MC events. 

The PDF for the B^ J/'?/'7r+7r~non-resonant signal is a parameterized M^+^- distribu- 
tion of the B^ — >■ J/tpTc'^TT^ phase space decay model. 

The AE distribution is used to determine the PDF for the background M^+^- shape. 
According to the MC study, background modes do not give a peak in the AE distribution. 
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Therefore, we estimate background by fitting the AE distribution with a first order polyno- 
mial (for background) and a Gaussian (for signals and K^; J/iI)tt^tt~ final state). In 
order to estimate the background as a function of Mn-+^~, we subdivide the candidate event 
sample in 0.2 GeV/c^ M^+^- bins from 0.25 GeV/c^ to 2.25 GeV/c^. Then a fit to the AE 
distribution in each M^+^- bin is performed. One of the major background sources here is 
the — s> J/ipK*'^, K*^ K^Ti~ decay, when a charged kaon is misidentified as a pion. Be- 
cause of the difference between the kaon and pion masses, this background mainly populates 
a region around AE ^ —0.1 GeV. By choosing the region of —0.05 GeV < AE < 0.2 GeV 
for our fits we avoid this contribution. The signal width and mean value are fixed according 
to signal MC expectations. Fig. |21 shows the results of the fits. 

We fit the obtained background with a fourth order polynomial as a function of M^^+t^-. 
Since this PDF is determined from data in a model independent way, we fix it in the fit. 

In order to find the probability of the B^ J/'^f2 and non-resonant B^ J/iIjtt^tt~ 
decays, the fit is carried out in three different cases: (i) PDFs for B^ J/ipp^, B^ J/ipKs 
and background are summed, (ii) the B^ contribution is additionally taken into 

account, and (iii) adding the non-resonant B^ — > J/ip-K^Tt" contribution to case (ii). The 
optimal likelihood values, £i, Cw and Cm for the cases (i), (ii) and (iii), respectively, are used 
to obtain a statistical significance: —2 In C\ = 42.4, —2 In = 33.7 and —2 In Cm = 33.7. By 
comparing cases (i) and (ii), a statistical significance of B^ — » J/4'f2 is found to be 2.9cr. As 
shown in Fig. El which is the fit result in case (ii), the M^+tt- spectrum is almost saturated 
by the Ks, and /2 contributions, and the non-resonant 5° — >■ J/t/'tt+tt" signal is found to 
be small. As a result of the fit in case (ii), the signal yields are obtained to be 139.8±17.0 for 
B^ J/ipp^ and 27.6±11.0 for B^ J /ipf2- To get a conservative upper limit for the non- 
resonant 5*^ — >■ J/ipTT^ir' contribution, the background PDF normalization is reduced by la 
and another fit to the M-^+j^- distribution is performed. The resulting yield is 15.1 it 20.1 
events. Taking systematic uncertainties listed in Table U and using the Feldman- Cousins 
approach [l3], an upper limit is obtained to be ;B(S° ^/'?/'(7r+7r^)non-rcs.) < 1 x 10^^ at 
90% C.L. 

For the B^ —>■ J/tpp^ and B^ —>■ J/'ipf2 modes, the following systematic uncertainties 
are taken into account: uncertainty of the track detection efficiency; background estima- 
tions; Breit-Wigner modeling, uncertainty of the relevant branching fractions; statistical 
uncertainties of the signal MC events; and uncertainty on the number of BB events. 

The tracking efficiency is estimated by MC and the control samples of partially recon- 
structed D* mesons, and its uncertainty is estimated to be 1.2% per track. The lepton 
identification efficiency is obtained from the ratio between the yields of single and double- 
tagged J/ip mesons, and its uncertainty is found to be 1.9% per track. An uncertainty 
of high-momentum charged pion identification is estimated from D*'^ D°(— K~tt^)7[~^ 
decays to be 1.7% per pion pair. The normalization of the background PDF is changed 
by ±lcr and a fit to M^^^ is performed again. Then the signal yield change is taken into 
account as a systematic error due to background estimation. An uncertainty due to L^s and 
R in the Breit-Wigner PDFs is estimated by the fit that is carried out changing L^s over 
its allower range. Two different values of the Blatt-Weisskopf barrier-factor radius R (0.5 
and 1.0 fm) are also compared. The maximum variation of the signal yield of each mode is 
assigned as a systematic error. Using a large MC sample of B'^ — *■ J/'^tt+tt" phase space 
decays, the detection efficiency is obtained as a function of M^+^- , and its product by the 
Breit-Wigner PDFs for and /2 mesons is integrated. The resulting values are compared 
to the detection efficiency obtained as an average over the whole kinematical region, and 
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the difference is assigned as a systematic error due to the signal MC modehng. 

All other investigated systematic uncertainties are small compared to those described 
above. We list all the estimated systematic uncertainties for each mode in Table lU 

TABLE I: Systematic uncertainties for each decay mode. Note that the uncertainty of background 
estimation for non-resonant J/i/jvr^vr^ is aheady taken into account to obtain a yield. 



Mode 




J/i^f2 


J/V'(7r+7r')non-rcs. 


Track reconstruction 


±4.8 % 


±4.8 % 


±4.8 % 


Lepton-ID 


±3.8 % 


±3.8 % 


±3.8 % 


7r±-ID 


±1.7 % 


±1.7 % 


±1.7 % 


Background estimation 


±5.1 % 


±15.4 % 




Breit-Wigner modeling 


±2.7 % 


±10.9 % 




Signal MC modeling 


±2.7 % 


±2.4 % 




MC statistics 


±0.7 % 


±0.7 % 


±0.7 % 


B{J/i) £+£-) 


±1.7 % 


±1.7 % 


±1.7 % 


B{f2 ^ vr+vr-) 




±2.9 % 




The number of BB 


±0.4 % 


±0.4 % 


±0.4 % 


Total 


±9.2 % 


±20.3 % 


±6.6 % 



Based on a large sample of BB pairs, we obtained the following preliminary branch- 
ing fractions for B decays resulting in the J/ipn^ir' final state: B{B^ J/ipp^) = 
(2.8±0.3(stat.)±0.3(syst.))xl0-5 and ^ = (9.8±3.9(stat.)±2.0(syst.)) x IQ-^ 

Since a statistical significance for the latter is 2.9a only, we also set an upper limit B{B^ 
J/i)f2) < 1.5 X 10-5 at the 90% C.L.; we also find B{B^ J/V'(vr+7r-) non-res.) < 1 x 10"^ 
at the 90% CL. The obtained value of B{B^ J jipp^) is consistent with that of BaBar 
and is more precise. 
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FIG. 2: The /S.E distribution for J/ijj'K^iT candidates in data in different M^+^- regions 

with the requirement 5.270 GeV/c^ < Mbc < 5.290 GeV/c^. In each plot the M^+^- range is 
0.2 GeV/c^: from (a) 0.25 GeV/c^ < M^+^- < 0.45 GeV/c^ to (j) 1.05 GeV/c^ < M^+^- < 
1.25 GeV/c^ in the alphabetic order. The combinatorial background is estimated by fitting this 
distribution with a Gaussian for a signal (the solid curve) and a first order polynomial for back- 
ground (the dashed line). 
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FIG. 3: The distribution of M^+^r- for — > J/ip7r~^iT^ candidates in data. The sohd hne shows 
the fitted results by taking the Kg, p^, f2 and background contributions into account. The dashed 
line is background only. 
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